function answer = polynom_i1(x, y, z)
 answer=x - 1.0295170327673958039493840033174*y - 0.029517032767395803949384003317391*x*z + 0.0006023884238244041622323265983141*x*((8*z)/3 - x*z + 224/3) + 0.0006023884238244041622323265983141*z*(10*x - 10*y);
 end
